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Abstract 



We investigate the structure of the constraints on three-point correlation functions 

tH- , emerging when conformal invariance is imposed in momentum space and in arbitrary 

^^ . space-time dimensions, presenting a derivation of their solutions for arbitrary scalar oper- 

^O ■ ators. We show that the differential equations generated by the requirement of symmetry 

'nJ" . under special conformal transformations coincide with those satisfied by generalized hy- 

^r^ . pergeometric functions (Appell's functions). Combined with the position space expression 

of this correlator, whose Fourier transform is given by a family of generalized Feynman 

(master) integrals, the method allows to derive the expression of such integrals in a com- 

kS ■ pletely independent way, bypassing the use of Mellin-Barnes techniques, which have been 

H ■ 
Cd ■ used in the past. The application of the special conformal constraints generates a new 

recursion relation for this family of integrals. 
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1 Introduction 

Conformal invariance plays an important role in constraining the structure of correlation func- 
tions of conformal field theories in any dimensions. It allows to fix the form of correlators - up to 
three-point functions - modulo a set of constants which are also given, once the field content of 
the underlying conformal field theory is selected |lll2]- The approach, which is largely followed 
in this case, is naturally formulated in position space, while the same conformal requirements, 
in momentum space, have been far less explored [3]. 

Conformal three-point functions have been intensively studied in the past, and a classification 
of their possible structures, in the presence of conformal anomalies, is available. Conformal 
anomalies emerge due to the inclusion of the energy momentum tensor in a certain correlator 
and, in some cases, find specific realizations in free field theories of scalars, vectors and fermions 
[Tl|2]. Typical correlators which have been studied are those involving the TT, TOO, TVV, and 
TTT, where T denotes the energy momentum tensor, V a vector and O a generic scalar operator 
of arbitrary dimension. 

The conformal constraints in position space, in this case, are combined with the Ward 
identities derived from the conservation of the energy momentum tensor and its tracelessness 
condition, valid at separate coordinate points, to fix the structure of each correlator. These 
solutions are obtained for generic conformal theories, with no reference to their Lagrangian 
realization which, in general, may not even exist. The solutions of the conformal constraints 
are then extended to include the contributions from the coincidence regions, where all the 
external points collapse to the same point. 

Free field theory realizations of these correlators (for fermions, scalars and vectors) allow to 
perform a direct test of these results both in position and in momentum space, at least in some 
important cases, such as the TVV or the TTT (this latter only for c? = 4) [3] , but obviously do 
not exhaust all possibilities. 

Recently interest in the momentum space form of conformal correlators has arisen in the 
context of the study of anomalous conformal Ward identities, massless poles and scalar degrees 
of freedom associated with the trace anomaly [H E], El [71 El El [10], and because of their possible 
role in determining the form of conformal invariance in the non- Gaussian features of the Cosmic 
Microwave Background [HI [12], or in inflation [131 [H]. The possibility of retrieving information 
on conformal correlators in momentum space seems to be related, in one way or another, to 
the previous knowledge of the same correlators in configuration space, where the conformal 
constraints are easier to implement and solve. One question that can be naturally raised is if we 
are able to bypass the study of conformal correlators in position space, by fixing their structure 
directly in momentum space and with no further input. This approach defines an independent 



path which, as we are going to show, can be successful in some specific cases. We will illustrate 
the direct construction of the solution bringing the example of the scalar three-point correlator. 
The analysis in momentum space should not be viewed though as an unnecessary complication. 
In fact, the solution in the same space, if found, is explicit in the momentum variables and 
can be immediately compared with the integral representation of the position space solution, 
given by a generalized Feynman integral. As a corollary of this approach, in the case of three- 
point functions, we are able to determine the complete structure of such an integral, which is 
characterized by three free parameters related to the scaling dimensions of the original scalar 
operators, in an entirely new way. It is therefore obvious that this approach allows to determine 
the explicit form of an entire family of master integrals. 

In the scalar case we are able to show that the conformal conditions are equivalent to partial 
differential equations (PDE's) of generalized hypergeometric type, solved by functions of two 
variables, x and y, which take the form of ratios of the external momenta. The general solution 
is expressed as a generic linear combination of four generalized hypergeometric functions of the 
same variables, or Appell's functions. Three out of the four constants of the linear combination 
can be fixed by the momentum symmetry. This allows to write down the general form of the 
scalar correlator in terms of a single multiplicative constant, which classifies all the possible 
conformal realizations of the scalar three-point function. 

In the final part of this work we go back to the analysis of the conformal master integrals, the 
Fourier transform of the scalar three-point correlators in position space. We show that the usual 
rules of integration by parts satisfied by these integrals are nothing else but the requirement 
of scale invariance. Specifically, dilatation symmetry relates the master integral J(z^i, z^2, ^'s), 
labelled by the powers of the Feynman propagators (z/i, 1^2, 1^3) - with ui + 1/2 + 1^3 = k, - to those 
of the first neighboring plane (/t — ?■ k + 1). On the other hand, special conformal constraints 
relate the integrals of second neighboring planes (/t — ?■ /t + 2). 

2 Conformal transformations 

In order to render our treatment self-contained, we present a brief review of the conformal 

transformations in rf > 2 dimensions which identify, in Minkowski space, the conformal group 

S0{2,d). 

These may be defined as the transformations x^ — )■ x' (x) that preserve the infinitesimal length 

up to a local factor 

dx^dx^ -)■ dx' dx'^ = VL{x)~'^ dx ^dx^ . (1) 



In the infinitesimal form, for d > 2, the conformal transformations are given by 

x'^{x) = x^ + a^ + ujyl^Xy + \x^ + h^x^ - 2x^h ■ x (2) 

with 

VL{x) = 1 — a{x) and a{x) = X — 2b ■ x . (3) 

The transformation in Eq.(|2]) is defined by translations (a^), rotations {Uf^^ = —Uu^), dilatations 
(A) and special conformal transformations (6^). The first two define the Poincare subgroup 
which leaves invariant the infinitesimal length {Q{x) = 1). 
If we also consider the inversion 

x,^x'^ = ^, nix)=x\ (4) 

^ x^ 

we can enlarge the conformal group to 0{2,d). Special conformal transformations can be 

realized by a translation preceded and followed by an inversion. 

Having specified the elements of the conformal group, we can define a quasi primary field 

O^x), where the index i runs over a representation of the group 0{l,d — 1) to which the field 

belongs, through the transformation property under a conformal transformation g 

0\x) 4 0'\x') = n{xfD]{g)0'{x) , (5) 

where r] is the scaling dimension of the field and DUg) denotes the representation of 0(1, d—1). 
In the infinitesimal form we have 

6gO'\x) = -{LgOy{x), with Lg = vd + r]a-^di,v,]W', (6) 

where the vector t>^ is the infinitesimal coordinate variation v^ = 6gX^ = x'^{x) —x^ and (S^^)* 
are the generators of 0(1, rf — 1) in the representation of the field O^ . The explicit form of the 
operator Lg can be obtained from Eq.([2]) and Eq.([3]) and is given by 

translations: Lg = a^d^ , 

rotations: Lg = — - [x^d^ - x^dy + S^^] , 

scale transformations : Lg = \[x ■ d + rj\ , 

special conformal transformations. : Lg = h^ {x^d^ — 2x^x ■ d — 2rjx^ — 2x^T,^''~\ . (7) 

Conformal invariant correlation functions of quasi primary fields can be defined by requiring 
that 

n 

Y^iOYi^i) ■ ■ ■ s,oir{xr) . . . ojr(x„)) = . (8) 

r=l 



In particular, the invariance under scale and special conformal transformations, in which we 
are mainly interested, reads as 



J2 (a^r ■ '9"'- + Vr) {0\' (a;i) . . . Olr{xr) . . . O'^M) = . 



r=l 



J2 {4d;^ - 22^./. a;, ■ d^^ - 2r]rXr^ - 2x,,(E/)};,) {Ol'{x,) . . . 0>(x,) . . . 0;"(x„)) = . 

(9) 



r=l 



The constraints provided by conformal invariance have been solved in coordinate space and for 
arbitrary space-time dimension. One can show, for instance, that the two and three-point func- 
tions are completely determined by conformal symmetry up to a small number of independent 
constants [IHl [1] . 

Exploiting the same constraints in momentum space is somewhat more involved. In the follow- 
ing we assume invariance under the Poincare group and we focus our attention on dilatations 
and special conformal transformations. 
For this purpose we define the Fourier transform of a n point correlation function as 

{2ny 6^''\pr + . . . +p„) {0\^{p,) . . . OJr(p„)) 

d'^xi . . . d^Xn {0\\xi) . . . C;"(x„))e*P^-"i+-+'P""", (10) 



where the correlation function in momentum space is understood to depend only on n — 1 
momenta, as the n-th one is removed using momentum conservation. 

The momentum space differential equations describing the invariance under dilatations and 
special conformal transformations are obtained Fourier-transforming Eq.Q. It is worth noting 
that some care must be taken, due to the appearance of derivatives on the delta function. As 
pointed out in [13j, these terms can be discarded and we are left with the two equations 

n-l 

" {OY{pi)...Oi:{Pr)...Ol-{pn))=0, 



n-l 



""1 / Q2 Q2 Q Q \ 

x(o'/(pi) ■ ■ ■ o';{vr) . . . o\:(p„)) = 0, (n) 

which define an arbitrary conformal invariant correlation function in d dimensions. Note that 
we are dealing with a first and a second order partial differential equations in n — 1 independent 
momenta. 



3 Two-point functions from momentum space and anoma- 
lies 

3.1 General solutions of the scale and special conformal identities 

We start exploring the implications of these constraints on two-point functions. In particular, 
the quasi primary fields taken into account are scalar (O), conserved vector (V^) and conserved 
and traceless (T^u) operators. 

For the two-point functions the differential equations in Eq. (iTTl) simplify considerably, being 
expressed in terms of just one independent momentum p, and take the form 

(-p^. ^+Vi + V2-d] G'^ip) = , 

where we have defined G^^{p) = {0\{p)02{—p))- The first of Eq. flT2l) dictates the scaling 

behavior of the correlation function, while special conformal invariance allows a non zero result 

only for equal scale dimensions of the two operators t]i = 772, as we know from the corresponding 

analysis in coordinate space. We start by illustrating this point. 

For the correlation function Gs{p) of two scalar quasi primary fields the invariance under the 

Poincare group obviously implies that Gs{p) = Gs{p'^), so that the derivatives with respect to 

the momentum p^ can be easily recast in terms of the variable p'^. 

The invariance under scale transformations implies that Gs{p'^) is a homogeneous function of 

degree a = ^{rji + ri2 — d) . At the same time, it is easy to show that the second equation in 

( 1T2|) can be satisfied only if r/i = ri2. Therefore conformal symmetry fixes the structure of the 

scalar two-point function up to an arbitrary overall constant C as 

Gsip') = {0,{p)02{-p)) = S,,,, C (p2)m-<^/2 . (13) 



If we redefine 



C - CS12 4,^_,/, ^^^^^ (14) 



in terms of the new integration constant csu, the two-point function reads as 
and after a Fourier transformation in coordinate space takes the familiar form 



Gsip') = K.. csu ::_„:'z:^^.' ip'r-"' , (is) 



(Oi(a:i)02(x2)) = 7.r. [Gsip")] = K,, 0312^^ , (16) 



where xu = xi — X2. The ratio of the two Gamma functions relating the two integration 
constants C and csu correctly reproduces the ultraviolet singular behavior of the correlation 
function and plays a role in the discussion of the origin of the scale anomaly. 

Now we turn to the vector case where we define Gy (p) = {V^' {p)V2 {—p)) ■ If the vector 
current is conserved, then the tensor structure of the two-point correlation function is entirely 
fixed by the transversality condition, S'^V^ = 0, as 

Gfip) = 7i'^^{p)Mp'), with ,r"/^(p) = ry-/^ - ^ (17) 

where fv is a function of the invariant square p^ whose form, as in the scalar case, is determined 
by the conformal constraints. Following the same reasonings discussed previously we find that 



with cy 12 being an arbitrary constant. Moreover we recall that a conserved current must have 
scale dimension rj = d — 1 for its divergence to be a quasi primary (scalar) field. To complete 
this short excursus, we present the solution of the conformal constraints for the two-point 
function built out of two energy momentum tensor operators which are symmetric, conserved 
and traceless 

T,,u = n^, 9%, = 0, T/ = 0. (19) 

Exploiting the conditions defined in Eq.( [T9|) we can unambiguously define the tensor structure 
of the correlation function G^ (]?) = 11'^ '^{p) /t(p^) with 



K'^'ip) = I 



t,^i^^p)t,P'^^p) + 7r"'^(p)7r^^(p) 



^ 7r"'^(p)7r^'^(p) , (20) 



d- 1 



and the scalar function frip"^) determined as usual, up to a multiplicative constant, by requiring 
the invariance under dilatations and special conformal transformations. We obtain 

Gr'^'ip) = Kv. CT.2 X^./^%~)'^'^ "^'"^^) ^^')"''^' • ^^^^ 

As for the conserved vector currents, also for the energy momentum tensor the scaling dimension 
is fixed and it is given hj t] = d. This particular value ensures that d'^T^^, is also a quasi primary 
(vector) field. 

These formulae agree with those in the literature p2j, and in particular those in Sec. 8 of Ref. 
[12] for the gravitational wave spectrum of the CMB. 



3.2 Divergences and anomalous breaking of scale identities 

The expressions obtained so far for the two-point functions in Eq.( fT5|) .( IT8|) and ( !2T|) . allow to 
discuss very easily the question of the divergences and of the corresponding violations that 
these induce in the scale identities. We can naturally see this noting that the Gamma function 
has simple poles for non positive integer arguments, which occur, in our case, when rj = d/2 + n 
with n = 0, 1, 2, . . .. 

Working in dimensional regularization, we can parametrize the divergence through an analytic 
continuation of the space-time dimension, d -^ d — 2e, and, then, expand the product r{d/2 — 
v) (p^)^^*^^^, which appears in every two-point function, in a Laurent series around d/2—r] = —n. 
We obtain 

r {d/2 - v) {j?r-"^ = ^ (-] + Hn + 1) + 0(e)) {p'r^^ , (22) 

where ip{z) is the logarithmic derivative of the Gamma function, and e takes into account the 
divergence of the two-point correlator for particular values of the scale dimension t] and of the 
space-time dimension d. 

The singular behavior described in Eq. (l22p is responsible for the anomalous violation of 
scale invariance [17], providing an extra contribution to the differential equation flT2l) obtained 
from the conformal symmetry constraints. Indeed, when rj = d/2 + n, employing dimensional 
regularization, the first of Eq. ( fT2l) becomes 



p^-^-n-e) G'^{p^) = 0, with r]^ = r]2 = v (23) 

which is the Euler equation for a function G^^ which behaves like (p^)"""*"^. Due to the appearance 
of a divergence in 1/e in the correlation function, Eq.( l23l) acquires an anomalous finite term in 
the limit e — )■ and we obtain 



d_ 

Qp2 



P" i-, - ^] GHp') = Gln,{p') . (24) 



where G^^^ [p^) corresponds to the singular contribution in the correlation function, which we 
have decomposed according to 



e 



As one can see from the r.h.s. of Eq.( p^ . the coefficient of the divergence, G^g^^gijP'), of the 
two-point function provides the source for its anomalous scaling. 



We illustrate the points discussed so far with some examples. Consider, for instance, the 
scalar correlator in Eq. flT5|) with scaling dimension rji = ri2 = rj = d/2. Due to the appearance 
of a pole in the Gamma function, the two-point correlator develops a divergence and becomes 



Gs{p ) = - C512 



where we have defined for convenience 



r {d/2) 



- + \ogp 



1 1 

- = -+1 
e e 



log(47r) 



(26) 



(27) 



with 7 being the Euler-Mascheroni constant. It is implicitly understood that the argument of 
the logarithm in Eq.( l26l) is made dimensionless, in dimensional regularization, by the insertion 
of a massive parameter. 

As one can easily verify, the scalar two-point function given in Eq. fl26|l satisfies the anomalous 
scaling equation fl2^ with a constant source term 



G 



S,sing 



ip' 



TT 



d/2 



C512 



(28) 



r {d/2) 

determined by the coefficient of the singularity. Note that the anomalous scaling behavior in 
Eq. flMl) is reproduced by the logarithmic contribution in Eq. (l26l) . 

Now we turn to the discussion of a correlation function with two vector currents. As already 
mentioned, the scaling dimension of the conserved vector operator is fixed at the value t] = d—1. 
In this case the divergences occur a.t d = 2n + 2 with n = 0, 1, . . . , so that, for d > 2, the first 
singularity appears at d = 4. Therefore the vector two-point function for d = 4 is 



Gf{p') 



^ 2 
Cyi2 — P 



- - 1 + logp 
e 



TC 



al3 



(P) 



(29) 



with e defined in Eq. (l27l) . As for the previous case, it is manifest that the two-point function 
in Eq. fl29|) satisfies the identity given in Eq. flMj) . with the logarithm accounting for the source 
of the anomalous scaling behavior. 

Finally, we illustrate the case of the correlation function built with two (symmetric, con- 
served and traceless) energy momentum tensors with scale dimension t] = d, which is slightly 
more involved, as we have to pay attention to the fact that 11'^ '^{p) itself depends on the space- 
time dimension d. The singularities are generated when d = 2n with n = 0,1, . . . , namely for 



even values of the space-time dimension. For instance, the two-point function in d 
to be given by 



TT 



G%^^^{p) = -CTu-:^{p') 



2\2 



'T 



192 



1 



+ logp^ 



<'"'{p)-7^^'''ip)^"'{p) 



4 is found 



(30) 



As we have already discussed previously, the appearance of the singularity in the correlation 
function develops an anomalous term in the scale identity. Correspondingly, being the energy 
momentum tensor related to the dilatation current, J{^ = x^T'^'^, it acquires an anomalous trace 
reflecting the violation of the scale symmetry. In this respect, the two-point function in Eq.( l30l) 
is characterized by a non vanishing trace 

V,. G-^^'^-'ip) = CTu ^ {p'r vr"^(p) , (31) 

generated by the last term in Eq. (!30|) which, on the other hand, arises from the explicit depen- 
dence of the n^ (p) tensor on the space-time dimension. The non-zero trace of Eq. (!3T]) is the 
signature of a conformal or trace anomaly, whose coefficients are known for free fields [HI [19] . 

4 Three-point functions for scalar operators 

In this section we turn to the momentum space analysis of conformal invariant three-point 
functions, by solving the constraints emerging from the invariance under the conformal group. 
We consider scalar quasi primary fields Oi with scale dimensions rji and define the three-point 
function 

Gl23(Pl,P2) = {Oiipi)02ip2)Osi-p, - P2)) . (32) 

The three-point correlator is a function of the two independent momenta pi and p2, from which 
one can construct three independent scalar quantities, namely pf, p\ and pi ■ p2. We trade the 
last invariant for p\ in order to manifest the symmetry properties of 6*123 under the exchange 
of any couple of operators. 

We observe that scale invariance, the first equation in Eq. fllll) . implies that 6*123 is a homo- 
geneous function of degree a = —d + |(?7i + ?72 + ^3)- Therefore it can be written in the 
form 

Gi23{pl,pl,pl) = {pl)-''-'-^''''-'^''''''^Hx,y) with x = 4, y = 4, (33) 



where we have introduced the dimensionless ratios x and y, which must not be confused with 
coordinate points. The dilatation equation only fixes the scaling behavior of the three-point 
correlator giving no further information on the dimensionless function $(x, ?/). 

The last equation of flTTl) . which describes the invariance under special conformal transfor- 
mations, is the most predictive one and, as we shall see, completely determines $(x,?/) up to a 
multiplicative constant. 
To show this, we start by rewriting Eq. fITT]) in a more useful form by introducing a change of 

10 



variables from (^1,^25^3) ^o (^^jZ/jPI)- The derivatives respect to the momentuni components 
are re-expressed in terms of derivatives of the momentum invariants and their ratios as 

d , 3 2 \ ^ r,f X y 9 



2{pl^,+P2^^)^-^ + — ((1 - x)pif, - xp2^,) t, 2(pi^ +p2/ 



ap'f opi^ PI ox Pi^y 

d d X 2, 

2(pi^ +^2^)7^:2 - '^iPl^^+P2^.)-2T^: + "2 ((^ - ^)P2m " VPif^) IT- (34) 



Similar but lengthier formulas hold for second derivatives. Also notice that the derivatives 
with respect to p| can be removed using the solution of the dilatation constraint in Eq. flS^ . 
Therefore we are left with a differential equation in the two dimensionless variables x and y. 
Due to the vector nature of the special conformal transformations, Eq. (fTT]) can be projected 
out on the two independent momenta pi and p2, obtaining a system of two coupled second 
order partial differential equations (PDE) for the function ^{x,y). After several non trivial 
manipulations, these can be recast in the simple form 



^{x,y) =0 



x{l- 


--)& 


7/2 9^ 
y dy^ 


-(a + /3 + l)4- 


-a/3 


y{i- 


y)dy2 




'^^VaZy + b' {a + P + ml 




\ 






-(« + /? + l)xj^- 


- a/3 



^{x,y) = 



(35) 



with the parameters a,/3,7,7' defined in terms of the scale dimensions of the three scalar 
operators as 

d r]i + rj2-r]z d 

a = - ^ , 7=2-^1 + 1' 

P = d , 7 = - - r/2 + 1 . (36) 

It is interesting to observe that the system of equations in (!35l) . coming from the invariance 
under special conformal transformations, is exactly the system of partial differential equations 
defining the hypergeometric Appell's function of two variables, F4(a,/3; 7,7'; a;, y), with coef- 
ficients given in Eq. (l36ll . The Appell's function F4 is defined as the double series (see, e.g., 
for thorough discussions of the hypergeometric functions and their properties) 



00 00 



F4(a,/3;7,7;x,y) = }_^}^ -- (37) 

j=o j=o w;« V / ;j ■ j- 

11 



where {a)i = r{a + i)/r{a) is the Pochhammer symbol. 

It is known that the system of partial differential equations fl35l) . besides the hypergeometric 

function introduced in Eg. ( 1371) . has three other independent solutions given by 

5*2(0, /3; 7, 7'; x,?/) = x^"'^F4(a-7 + l,/3-7 + l;2-7,7';x,y), 
S3{a,l3;i,7';x,y) = y^"^' F4(a - 7' + 1, /3 - 7' + 1; 7, 2 - 7'; x,y) , 
S,{a,l3;^,^';x,y) = x^"^ y^-^' F4(a - 7 - Y + 2, /? - 7 - V + 2; 2 - 7, 2 - 7'; x,y) . 

(38) 

Therefore the function $(a;, y), solution of (|35|1 . is a linear combination of the four independent 
hypergeometric functions, i.e. 

Gi23{pIpI,pI) = {plr'-''-^'''-''''-'''''^Hx,y) 

4 

i=l 

where we have denoted with Si the Appell's function F4 given in Eq. (l37p . while the parameters 
a,/3,7,7' are defined in Eq. fl36|) . The Ci{rji,rj2,rj^) appearing in the linear combination, are the 
arbitrary coefficients which may depend on the scale dimensions rji of the quasi primary fields 
and on the space-time dimension d. 

The coefficients Ci{rii,ri2,r]3) can be determined, up to an overall multiplicative constant, by 
exploiting the symmetry of the correlation function under the interchange of two of the three 
scalar operators present in the correlator, which consists of the simultaneous exchange of mo- 
menta and scale dimensions {pj^rji) <-)■ {p'^,rij). 

Consider, for instance, the invariance of the three-point function under the exchange 02(^2) ^ 
Os{—pi —P2), which is achieved by {P2,V2) ^ {phV3)- Then Eq. fl39|) becomes 



~ X 1 

.Ci{rn,V3,V2)^i\'' "'''' ""' 

i=l 

where 



Gi32{pIpIpI) = (p^)-'^+^(''^+^'^+''^)$^Q(r;i,7^3,^2)5.(«,A7,7;-,-), (40) 

^^ y y 



~ , ^ d Vi + Va ~ V2 ~ , . d 

a = a{r]2 ^ m) = , 7 = 7(772 ^ t]^) = -- t]i + 1 = -f , 

P = I3{V2 ^ V3) = d- !^1±^±:^ = ;3, 7' = 7'(r/2 ^m) = ^-V3 + l- (41) 

Note that the hypergeometric functions are now evaluated in x/y and 1/y. To reintroduce the 
dependence from x and y, in order to exploit more easily the symmetry relation 

Gi23{pIpI,pI) = GMpIpIpD , (42) 

12 



we make use of the transformation property of F4 [20] 

F4(«,/3;7,y;x,y) = ^p^)^itlA(^y)-- F,{a,a - ^' + l-^,a - f3 + I--,-) 

r(Y-a)r(/3) y y 

+ S^?^^%^(-^)"' ^4(/3 -7' + l,/3; 7, /3 -« + !;-,-)• (43) 
r(7 -P)r(a) y y 

After some algebraic manipulations, and repeating the procedure described so far for the other 
operator interchanges, the ratios between the coefficients q take the simplified form 

c^{vi,V2,v3) r(f-r^,)r(-f + f + f)r(f-f + f-f) ' 

C2{vuv2,vs) ^ r {r^,- l)r if -f + f)r il + f -f -f) 
c,{vi,v2,v3) r (f - r/2) r (f + f - f ) r (-f + f + f + f ) ' 

ci{vuV2,v3) ^ r (r^, - f ) r (r^, - f ) r (d - f - f - f ) r (f - f - f + f ) 
c.{vuV2,V3) r{l-vi)r{l-v2)T{f + f-f)T{-l + f + f + f) 



(44) 



and define ^123(^1,^2)^3) ^P ^^ ^ multiplicative arbitrary constant C123 = £123(771, 772, '73)- This 
depends on the space-time dimension d, on the scale dimensions rji of the quasi primary fields 
and on their normalization. 

The conformal invariant correlation function of three scalar quasi primary fields with arbitrary 
scale dimensions is then given by 

2 2 2>, C123 7r'^4^"^('''+''2+''^)(pi)"'^+^(''^+'''+''') 

123\Pi,P2jP3) — -p /2j_ 1 22 _ 23\ T (m. _ 22 I 22. ^ r ( — Hi. _l 22 i V3\ -p ( _d i 21 _l 22 i 22. ^ 
-^l2"'"2 2/'-^l2 2"i"2Ml 2"'"2"i"2Ml 2"i"2"'"2"i"2/' 

r(„-0r(,«-0r(.-|-|-|)r(^-|-|.| 

/d m + m-m . vi + m + m d d 

I 2 2 ' 2 ^2 ~^^ '2 ~^ '^'^ 

ir/^*^ Ap/^ d\rfm V2,V3\j.fd Vi V2 % 

.r(.-0r(^-.)r(-|.|.|)r(^-|.|-| 

.rj2-i t:. fd VI+V3-V2 V2 + V3-Vl^d , ^ d 



xy''" 2F4 ( ^ , ■,--7]i + l,-- + 7]2 + l;x,y 

/d \j.fd \^fVi,V2 V3\ J. f d vi V2 ri3 

' ' (45) 
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The convenient normalization employed in Eq. fHSj) for the three-point function reproduces, 
through the operator product expansion, as we are going to show next, the normalization of 
the two-point functions which we have chosen in Eq.( TT5|) . 

As we shall identify the three-point correlator discussed in this section with specific Feynman 
amplitudes, this will fix the arbitrary constant C123 using some information coming from the 
same operator product expansion analysis. This topic will be presented in section |5l Indeed, the 
solution of the momentum space version of the conformal constraints provides an alternative 
computational tool for correlation functions with conformal symmetry. 

It is worth to emphasize the connection between the invariance under special conformal 
transformations and appearance of the Appell's functions. Indeed we have shown how the 
constraints provided by the conformal group translate, in momentum space, in the well-known 
system of partial differential equations defining the hypergeometric series F4. We have analyzed 
this connection in the case of a conformally invariant three-point function built with scalar 
operators in some detail. A similar correspondence should also hold for more complicated 
vector and tensor correlators. 

4.1 The Operator Product Expansion analysis 

In this section we show the consistency of our result with the operator product expansion (OPE) 
in conformal field theories in which the structure of the Wilson's coefficients is entirely fixed 
by the scaling dimensions of the two operators. 
Considering, for instance, the coincidence limit in the scalar case, one has 

0,{x^)0,{x2)r^y2-—^ -Okix2) for x,^X2, (46) 

where x^ = Xi — 0:2. It is worth noting that the coefficients Cijk are the same structure constants 

appearing in the three-point functions. 

For the correlation function of three scalar operators the OPE implies the singular behavior 

{0,{x,)0,{x,)Os{x,)) ^'^^^ J^l^l iO,{x,)0,{x,)) , (47) 

(^23)^ 

with analogous formulae for the other coincidence limits. For the sake of simplicity, we chose a 
diagonal basis of quasi primary operators normalized as 

(0.(xi)0,(x2)) = -|^. (48) 
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The momentuiii space version of the OPE in Eq. fH7|) reads 

d/2 Y(- — Hl±IMzin) p 

where the scalar two-point function is normahzed as in Eq.f llSp with C512 = 1. In the previous 
equation the symbol ~ stands for the momentum space counterpart of the short distance limit 
3^3 -^ ^2 which is achieved by the pl,P2 -> C)0 limit with P2/P3 ~^ ^■ 

The result for the scalar three-point function given in Eq.f l45p is indeed in agreement, as ex- 
pected, with the OPE analysis. This can be shown from Eq.f l45l) by a suitable expansion 
of the corresponding Appell's functions. In particular, in order to reproduce the momentum 
space singular behavior of Eq. (H9!) , we need the hypergeometric leading expansion in the limit 
X = pI/pI — )■ and y = pi /pi — )■ 1, which reads as 



rcvircv -a- b) 

F4(a,/3;7,y;x,i/)^ V ^\ ^ ^ for x^O,y^l. (50) 

r(7' - a)r(7' - fj) 

In the previous equation we have retained only the terms with the correct power-law scaling in 
the pI variable, as dictated by the OPE analysis. In this case these contributions come from 
the terms of Eq.( H5|) which are proportional to the 5*2 and 5*4 solutions defined in Eq.(l38l). 
Analogously, in the limit pl,pl — )■ 00, with pf/pl — )■ 1, which is described in coordinate space 
by X3 — !■ Xi, the leading behavior is extracted from 5*3 and ^4. 

The remaining coincidence limit Xi — ?■ X2, corresponding to pl,pl -^ 00 with p\/p\ — )■ 1, is more 
subtle due to the apparent asymmetry in the momentum invariants Pi,P2,P3 of the three-point 
scalar correlator, as given in Eq. (l45l) . In this case both x and y grow to infinity while their 
ratio x/y — t- 1. Therefore it is necessary to apply the transformation defined in Eq. fj43|) to each 
hypergeometric function appearing in Eq. fj45|) . This can be viewed as an analytic continuation 
outside the domain of convergence |a/x| + \y/y\ < 1, where the Appell's function is strictly 
defined as a double series. The hypergeometric functions are then expanded according to 

p. . , ^ . .-a r(7)r(y)r(/3-a)r(7 + y-2a-i) 

F4(a,/^;7,7;x,?/) ~ {-y) 



r(/3)r(7 - a)r(Y - a)r(7 + y - a - 1) 



_g r(7)r(y)r(a-/3)r(7 + y- 2/3-1) 



r(a)r(y - /3)r(7 - /3)r(7 + y - /? - 1) ' 

X 

for x,y ^ 00 , yl. (51) 

This completes the analysis of the OPE on the three-point scalar function in the three different 
coincidence limits. 
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5 Feynman integral representation of the momentum 
space solution 

We have seen in the previous sections that we can fix the exphcit structure of the generic 
three-point scalar correlator in momentum space by solving the conformal constraints, which 
are mapped to a system of two hypergeometric differential equations of two variables. These 
variables take the form of two ratios of the external momenta. In particular we find that in any 
d dimensional conformal field theory the solution of this system of PDE's is characterized by a 
single integration constant which depends on the specific conformal realization, as expected. 

In this section we want to point out the relationship between the scalar three-point functions 
studied so far and a certain class of Feynman master integrals. These can be obtained by a 
Fourier transformation of the corresponding solution of the conformal constraints in coordinate 
space, which is well known to be 

{0,{x,)0,{x2)0s{xs)) = ^- ^ ^ ^ rr- 7 • (52) 

Transforming to momentum space, we find an integral representation, which necessarily has 
to coincide, up to an unconstrained overall constant, with the explicit solution found in the 
previous section, and reads as 

f d'^l 1 

with external momenta pi, p2 and ps constrained by momentum conservation pi + P2 + P3 = 
and the scale dimensions rji related to the indices z/j as 

ril = d-V2~V^, V2 = d - 1^1 - 1^3 , 7^3 = d - z/i - z/2 . (54) 



This expression describes a family of master integrals which has been studied in [23l UM ; whose 
explicit relation with Eq.f lS^ is given by 



(i°'pi (i°'p2 d^ps 



{2'Ky6^^\p^+P2+P3)J{iyi, 1^2,^3] 



{2nY {2nY {27tY 

1 T{d/2 - u^)T{d/2 - U2)T{d/2 - u^) 



e 



-tpi-Xi—tp2-X2—ip3-X3 



4u,+U2+U3^3d/2 r(z/i)r(z/2)r(z/3) (x22)'^/2-'^3(a;2^)<i/2-^i(a;2Jd/2-^2 ' 

(55) 

The integral in Eq.( l53l) satisfies the system of PDE's ( l35l) . Therefore, it can be expressed 
in terms of the general solution given in Eq. (l45|) which involves a linear combination of four 
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Appell's functions, with the relative coefficients fixed by the symmetry conditions on the de- 
pendence from the external momenta. Then Eg. (145 p identifies J(z/i, 1^2,^3) except for an overall 
constant C123 which we are now going to determine. This task can be accomplished, for instance, 
by exploiting some boundary conditions. 

As for the OPE analysis discussed in the previous section, we may consider the large momentum 
limit in which the three-point integral collapses into a two-point function topology. Taking, for 
instance, the P2,pl -^ 00 limit with p^/pl — ?■ 1 we have 

A-^.-^. "=) ~ (^ / (g. (,.).((/+ p.).). = (^(£j^ eC'- ''3) (Pif^-"- . (56) 
where 

Eq.f l56p must be compared with the same limit taken on the explicit solution in Eg. (145 p . where 
the scale dimensions rji are replaced by z/j through Eq. fl5^ . This completely determines the 
multiplicative constant C123 and the correct normalization of the three-point master integral, 
which is obtained by choosing 

^ z^-'^ T{d/2 - u^)T{d/2 - z/2)r(^/2 - z/3) . . 

C123 4,,+,,+,3^3d/2 r(z/i)r(z/2)r(z/3) ' ^ ' 

Therefore the scalar master integral is given by 

Jii^i, 1^2, 1^3) = G 123i.pl, pI^pI) (59) 

with scaling dimensions defined in Eq. (l54p and the coefficient C123 in Eq.f l58|) . Notice that 
this method allows us to bypass completely the Mellin-Barnes techniques which has been used 
previously in the analysis of the same integral. 

5.1 Recurrence relations from conformal invar iance 

Having established the conformal invariance of the generalized three-point master integral 
</(^i, ^2, 1^3)1 we can study the implications of the conformal constraints on the integral repre- 
sentation of Eq.( [53|) . These are automatically satisfied by the explicit solution given in Eq.( H5|) . 
but once that they are applied on J(z/i, 1/2, z/3), generate recursion relations among the indices 
of this family of integrals. Specifically, they relate integrals with ui + 1/2 + 1^3 = k to those with 
i^i + 1^2 + ^3 = 1^ + ^ and z/i + z/2 + z/3 = K + 2. For instance, differentiating Eq.(j53ll under the 
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integration sign according to the first of Eq.f llip . wliicli is tlie condition of scale invariance, we 
easily obtain the recursion relation 



Z/2P1 J(z/i, V2 + 1, y^) + ^iPI Ji^l + 1, ^2, v-i) 



{vi + 1^2 + 2 z/3 - d) J(z/i, z/2, z/a) 



+ V2 J{Vl, V2 + 1, Z/3 - 1) + Vi J{Vi + 1, z/2, Z/3 - 1) 



(60) 



together with the corresponding symmetric relations obtained interchanging {p\^ vi) ^ (pg, z/3) 
or (p|, 1^2) ^ ij>li ^^)- These equations link scalar integrals on two contiguous planes, as men- 
tioned above. The recurrence relations obtained from scale invariance exactly correspond to 
those presented in [24j and following from the usual integration-by-parts technique, which in 
this case is derived from the divergence theorem in dimensional regularization 



dH d 



I, 



{2'ky di, WYK{i+pify-{{i-P2?y^ 



0. 



(61) 



We can easily show the equivalence between Eq. ( 16111 and the first of Eq.l jTTjl which is the 
constraint of scale invariance. In fact, the scale transformation acts on J(z/i, z/2, z/3) in the form 



d-2 (z/i + z/2 + z/3) - pi ■ p2 ■ 7^— 

opi ap2 



dH 



1 



{py^ {{I + PifY^ {{I - P2fy^ 



0. (62) 



Now we just invoke Euler's theorem on homogeneous functions on the integrand, which is of 
degree —2 (z/i + z/2 + z/3) in the momenta pi, p2 and I and obtain the relation 



d d ^ d 

Pi • ^— + P2 ■ ^— + / ■ ^ 

Opi Op2 Ol 



1 



{py'{{i+Piyy^{{i-P2?y- 

-2(z/i + z/2 + z/3) 



{i^y'{{i+Piyy^{{i-P2?y^ 



■ (63) 



At this point, if we combine Eqs. (l62p and (l63l) and rewrite (i as ^ ■ Z, we easily obtain the 



equivalence with Eq. fl6T|) . 

Other recursive relations can be found requiring Eq. fl53|) to satisfy the constraint of special 
conformal invariance which, from the second equation in Eq. fITT]) . takes the form 



Pi/ 



92 



dpi ■ dpi 



-2pi 



92 



''A^M 



dp'^dpi y 



2 (z/2 + z/3)^ + (1 o 2) [^ J(z/i, z/2, z/3) = . (64) 



This is a vector condition which involves some tensor integrals of the same J(z/i, z/2, t's) family. 
Differentiating the integral J(z/i, z/2, z/3) as in Eq. flMj) and performing some standard manipula- 
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tions one arrives at the implicit formula 

^2PltM (1 + 2^2 + 1^3 -d/2) J{vi,V2 + 1,1^3) + (z^2 + l) (^(^^1,2^2 + 2, I/3 - 1) -p? J(j^l,Z^2 + 2, Z^s)) 

+ ^^lP2^l (1 + V1 + V2. -d/2) J{vi + 1,1^2,1^3) + (1^1 + 1) (-''(1^1 + 2,z^2,i^3 - 1) -p^Ji'^i + 2,2^2,^^3)) 

+ U2 (z^3 - l)Jf,{iyi,U2 + 1,1^3) + (z^2 + 1) (^/i(l^l,J^2 + 2, 1^3 - 1) - Pi J^(z^l, Z^2 + 2,1^3)) 

(z^3 - 1) Jf,{iyi + 1, 1/2, J^s) + (i^i + 1) (-''/^(i^i + 2, 1/2, J/3-1) - P2 -^^(^1 + 2, z^2, J^s)) 



-Z/1 

where the rank-1 tensor integral is defined as 
d'^l I,. 



0, (65) 



J^,{l^l,J^2,J^3) = I 



Ci{ui, Z/2, Z/3) Plf,- 02(1^1, Z/2, 1^3) P2f, , 

(66) 



(27r)'^ (/2).3((/+p,)2).2((/_p2)2).i 

with the coefficients given by 

C^{U„U2,U3) = (^2_p2_p\)2_4p2^2 {(Pl+P2-pD^(^l-l,^2,^3) 

-2p2 ^(Z^l, Z^2 - 1, 1^3) + (-Pl + P2 + Ps) Ji^l^ ^2, Z^3 - 1) + P2 (Pl -P2+ pI) Ji.^1^ ^2, ^^s) \ 



C2(z/i,Z/2,Z/3) 



1 



(P? + P2 - pD Ji.^1^ ^^2 - 1, ^^3) 



(p2_p2_p2)2_4p2p2 

-2pl J(z/i - 1, 1/2, Z/3) + (pi - P2 + P3) Ji.^li ^2, Z/3 - 1) + P? (-Pi + P2 + P3) J{^1, ^2, 1^3) > ■ 

(67) 



Using the momentum expansion of the tensor integral defined above, we extract from Eq.fl65 
the relation 

z/2 (1/3 - 1) Ci(l/i, Z/2 + 1, 1/3) + Z/2 (1/2 + 1) (Ci(l/i, Z/2 + 2, Z/3 - 1) - p^ Ci(l/i, Z/2 + 2, Z/3)) 

-z/i (z/3 - l)Ci(z/i + 1, z/2, 1/3) -z/i (z/i + 1) (Ci(z/i + 2, z/2, 1/3 - 1) -p2C'i(z/i + 2, 1/2, z/3)) 

(z/2 + 1) (J(Z/1, Z/2 + 2, 1/3 - 1) - pj J(Z/1, Z/2 + 2, 1/3)) + (1 + z/2 + z/3 - (i/2) J(Z/1, z/2 + 1, I/3) 



+ Z/2 



= 0, 

(68) 



together with the corresponding symmetric equation obtained interchanging {pi, z/i) ^ {p\,V2)- 
This result allows to express integrals in the plane z/i + z/2 + z/3 = k + 2 in terms of those in the 
two lower ones. In fact, introducing in Eq. (l68p and in its symmetric one the explicit expressions 
for Ci and C2 we get 

1 



J(z/i + 2,z/2,z/3) 



yi{vi + l){pl+pl-pl)plpl 



J]] C(a,6,c)^(z^i + a, z/2 + 6, z/3 + c), (69) 



(a,h,c) 
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where the coefficients C(a,b,c) are given by 



C{o,o,o) = (^^3 - 1) ( («^i + 1^2) Pi-i^2P3J , 

C{i,-i,o) = J^i {J^3 - 1) {pI - pD , 
C{_i,i,o) = -1^2(1^3- l)p?, 

C{o,i -1) = 1^2 {^2 + 1) Pi - (2 + z/2 - z/s) P3 
C(i,o,-i) = 1^1 {pli^i + 1) - Psl'^s - 1)) , 

C{2,-i,-i) = -«^i («^i + 1) {pI-pD , 

C(-l,2-l) = -//2(«^2 + l)Pl, 

C{2,0-2) = -Z/l(z/l + 1)^3, 

C{0,2-2) = V2{V2 + '^)pI, 
'd 



C(1,0,0) = ^1 



{pi 



2n2 



2 



rf 



P3 {P2 - Pa) o - ^1 - ^3 - 1 



+ Pf( (l--)p^+P^(2z/i + z/3 + 3-rf) 



C(0,1,0) = ^2Pi 

C 



[l-^){pl-pl) + i^-2-u2)pl 



(0,2,-1) = -J^2 {J^2 + 1) Pi P3 , 

C(-l,2,0) = t^2 {V2 + 1) (P?)^ ) 

C{2,o,-i) = z^i (z^i + 1) P3 (Pi + 2pl-pl), 

C(2,-l,0) = «^1 («^1 + 1) P2 (Pl - P3) • 



(70) 



Analogous results hold for J{vi^V2 + 2, z/3) and J(z^i, z>'2, ^^3 + 2) if we just make the usual 
exchanges {pi, ui) ^ {p\, U2) and {pi, Ui) ^ {p\, u^) both in the integrals and in the coefficients 



C, 



{a,b,c)- 



6 Conclusions 

We have shown that the solution of the conformal constraints for a scalar three-point function 
can be obtained directly in momentum space by solving the differential equations following 
from them. This has been possible having shown that these constraints take the form of a 
system of two PDE's of generalized hypergeometric type. The solution is expressed as a linear 
combination of four independent Appell's functions. The use of the momentum symmetries of 
the correlator allows to leave free a single multiplicative integration constant to parameterize 
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the general solution for any conformal field theory. If this solution is compared with the posi- 
tion space counterpart and its Fourier representation, which is given by a family of Feynman 
master integrals, we obtain the explicit expression of the same integrals in terms of special func- 
tions. Our solution coincides with the one found by Boos and Davydychev using Mellin-Barnes 
techniques, which in our case are completely bypassed. Having established this link, we have 
shown that by applying special conformal constraints on the master integral representation one 
obtains new recursion relations. 

The momentum space approach discussed in this work can be used to treat more complicated 
correlators. For instance, this method can be employed in the analysis of three-point func- 
tions involving the vector and the energy momentum tensor operators, like VVV, TOO, TVV 
and TTT, as well as higher order ones, such as the scalar four-point function, whose general 
structure has been known for a long time [I5]. Nevertheless, such a treatment is much more 
complicated, in the former case due to the tensor nature of the correlators, which implies a much 
more involved set of constraints, in the latter because of the increasing number of independent 
variables in the partial differential equations. This is left for future work. 

References 



[1] 

[2] 
[3] 

[4] 
[5] 
[6] 

[7] 



H. Osborn and A. C. Petkou, Ann. Phys. 231, 311 (1994), arXiv:hep-th/9307010 



J. Erdmenger and H. Osborn, Nucl.Phys. B483, 431 (1997), arXiv:hep-th/0103237 



C. Coriano, L. Delle Rose, E. Mottola, and M. Serino, JHEP 1208, 147 (2012), 
larXiv:1203.1339. 

M. Giannotti and E. Mottola, Phys. Rev. D79, 045014 f2009). larXiv:0812.0351[ 

E. Mottola, Acta Physica Polonica B41, 2031 ^2010). [arXiv:1008.5006[ 

R. Armilhs, C. Coriano, and L. Delle Rose, Phys. Rev. D81, 085001 (2010), 
larXiv:0910.33811 

R. Armilhs, C. Coriano, and L. Delle Rose, Phys. Lett. B682, 322 f 2009). ErXTv: 0909. 45221 

R. Armillis, C. Coriano, and L. Delle Rose, Phys.Rev. D82, 064023 (2010), 
larXiv: 1005.41731 

[9] C. Coriano, L. Delle Rose, A. Quintavalle, and M. Serino, Phys.Lett. B700, 29 (2011), 
larXivillOl. 16241 



21 



[10] C. Coriano, L. Delle Rose, and M. Serine, Phys.Rev. D83, 125028 f2011). la7XiV:1102.4558[ 
[11] I. Antoniadis, P. O. Mazur, and E. Mottola, Phys. Rev. Lett. D79, 14 (1997), 



astro-ph/9611208. 



[12] I. Antoniadis, P. O. Mazur, and E. Mottola, JCAP 09, 024084 r2012). iFxrv:1103.4164[ 

[13] J. M. Maldacena and G. L. Pimentel, JHEP 1109, 045 ^20111. ETXTv: 1104.28461 

[14] A. Kehagias and A. Riotto, Nucl.Phys. B864, 492 (2012), arXiv: 1205. 15231 

[15] A. M. Polyakov, JETP Lett. 12, 381 (1970). 

[16] A. Cappelli, D. Friedan, and J. Latorre, Nucl.Phys. B352, 616 (1991). 

[17] L. S. Brown and J. C. Collins, Ann. Phys. 130, 215 (1980). 

[18] M. J. Duff, Nucl. Phys. B125, 334 (1977). 

[19] P. R. Anderson, C. Molina-Paris, and E. Mottola, Phys. Rev. D67, 024026 (2003), 0209075. 

[20] P. Appell and J. Kampe de Feriet, Fonctions hypergeometriques et hyperspheriques 
(Gauthier-Villars, 1926). 

[21] W. N. Bailey, Generalized hypergeometric series (Cambridge University Press, 1935). 

[22] L. J. Slater, Generalized hypergeometric functions (Cambridge University Press, 1966). 

[23] E. Boos and A. I. Davydychev, Moscow Univ.Phys.Bull. 42N3, 6 (1987). 

[24] A. I. Davydychev, J.Phys.A A25, 5587 (1992). 



22 



